FERMILAB - TM - 2529 - E 



The DO Run lib Luminosity Measurement 

B.C.K. Casey, ^ M. Corcoran/ K. DeVaughan/ Y. Enari/ E. Gallas*/ I. Katsanos/ J. Linnemann/ 
J. Orduna/ R. Partridget/ M. Prewitt/ H. Schellman/ G.R. Snow/ and M. Verzocchi^ 

'LPNHE, Universites Pans VI and VII, CNRS/IN2P3, Pans, France 
^Fermi National Accelerator Laboratory, Batavia, IlUnois 60510, USA 
'''Northwestern University, Evanston, Illinois 60208, USA 
^Michigan State University, East Lansing, Michigan 48824, USA 
^University of Nebraska, Lincoln, Nebraska 68588, USA 
^Brown University, Providence, Rhode Island 02912, USA 
''Rice University, Houston, Texas 77005, USA 
(Dated: April 02, 2012) 

An assessment of the recorded integrated luminosity is presented for data collected with the 
DO detector at the Fermilab Tevatron Collider from June 2006 to September 2011 (Run lib). In 
addition, a measurement of the effective cross section for inelastic interactions, also referred to 
as the luminosity constant, is reported. This measurement incorporates new features that lead 
to a substantial improvement in the precision of the result. A luminosity constant of ctlm ~ 
48.3±1.9±0.6 mb is obtained, where the first uncertainty is due to the accuracy of the inelastic cross 
section used by both CDF and DO, and the second uncertainty is due to DO sources. The recorded 
luminosity for the highest Et jet trigger is £i.cc = 9.2 ± 0.4 fb~^, with a relative uncertainty of 
4.3%. 



PACS numbers: 13.85.Lg 

I. INTRODUCTION 

An essential ingredient in cross section measurements 
is the integrated luminosity, £, used to normalize the 
data sample. At DO, the instantaneous luminosity, L, 
is derived from hit rates produced from inelastic proton 
- antiproton collisions registered in a dedicated detector 
system. Measured hit rates are converted to luminosity 
using a normalization procedure based on the total in- 
elastic cross section, and the geometric acceptance and 
efficiency of the dedicated detector system for register- 
ing inelastic events. The measurement of the effective 
cross section for inelastic interactions, and the assessed 
recorded integrated luminosity for data collected with the 
DO detector at the Fermilab Tevatron Collider from June 
2006 to September 2011 are reported. Luminosity L var- 
ied during that period in the range (5 - 420) /ib^^s^^ 
(equivalent to (5 - 420)-1030 cm^^g-i). 

In this note, a short description of the detector used 
for the instantaneous luminosity measurement and of the 
data samples used for this study is followed by a review 
of the luminosity measurement technique. The following 
sections describe the backgrounds that affect the lumi- 
nosity measurement, the calculation of the detector ac- 
ceptance, and the calculation of the luminosity constant 
and its uncertainty. In the Appendices, the luminosity 
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measurement technique and the background removal are 
described in more detail. 

A. Luminosity Monitor Detector 

The Luminosity Monitor (LM) [l|, Q consists of two 
arrays of scintillation counters mounted on the DO end- 
cap calorimeter cryostats as indicated in Fig. [1] In the 
description of the DO detector a right-handed coordinate 
system is used. The z-axis is along the proton beam 
direction. The angles (j) and are the azimuthal and 
polar angles, respectively. The r coordinate denotes the 
perpendicular distance from the z axis. 

From the perspective of the proton beam, the upstream 
LM array is called the "north" LM and the downstream 
array is called the "south" LM. Each array has 24 wedge- 
shaped scintillation counters with fine-mesh photomulti- 
plier tube (PMT) readout. The PMT signals are ampli- 
fied on the detector and are carried on low-loss cables to 
the LM VME electronics where the charge and the tim- 
ing of PMT signals are measured. Coverage is provided 
over the pseudorapidity interval of 2.7 < \r]\ < 4.4, where 
ri = -ln[tan(6l/2)]. 

The LM electronics identify in-time hits that are within 
±6.4 ns of the nominal time-of-flight from the center of 
the DO detector to the LM. This window is about three 
times the width of the time distribution for in-time hits. 
Halo particles typically produce hits that are ^9 ns early 
in one of the detectors. A luminosity coincidence is iden- 
tified when there is at least one in-time hit in both the 
north and the south LM detector arrays. Since beam 
crossings with many early hits from beam halo interac- 
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FIG. 1: The Luminosity Monitor layout. In (b) the sohd dots 
represent the location of the PMTs. 



tions can lead to luminosity measurement errors, a "halo 
veto" is applied when there are six or more early hits in 
one or both detector arrays. The beam crossings that do 
not trigger the halo veto are called "live crossings"^. 



B. Data Samples 

The data sample that DO recorded during Run II of 
the Fermilab Tevatron Collider is split in two periods: 
(i) data collected between April 2002 and February 2006 
(Run Ila), and (ii) data collected from June 2006 to 
September 2011 (Run lib). A major difference between 
the two periods is the addition of an inner silicon layer Q 
to the DO Sihcon Microstrip Tracker [i] (SMT) during the 
2006 shutdown. Other differences between the two peri- 
ods include removal of a forward silicon disk on each end 
of the SMT and introduction of a new beryllium beam 
pipe with a flange near the LM. The readout system of 
the LM detector was upgraded between Runs lla and lib 
to reduce the electronic noise jH . 

The data from the LM detector information includes 
measurements of the arrival time and pulse height in- 
formation for each of the 48 LM counters. In addition, 
the LM electronics allow the accumulation of histograms 
of quantities calculated by the LM electronics for cali- 
bration and monitoring purposes. These histograms are 
accumulated at the beam crossing rate with no deadtime. 

The Fermilab Tevatron Collider has 1113 possible radio 
frequency (RF) buckets. The minimal spacing between 
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FIG. 2: 2D multiplicity distribution, after background sub- 
traction, collected with the histogramming feature of the LM 
electronics. 



RF buckets where particles can be placed is one "tick" 
and corresponds to a gap of 132 ns. One turn of the Teva- 
tron consists of 159 ticks, 36 of which generally contain 
beam. The ticks that actually contain particles are called 
"beam bunches" , and the collision of proton and anti- 
proton bunches is called a "beam crossing" or "bunch 
crossing" . The beam bunches are arranged in 3 evenly 
spaced "bunch trains" , separated by a 2.5 /xs abort gap, 
and within each bunch train there are 12 beam bunches, 
each separated by 396 ns. Ticks that do not contain beam 
are referred to as "empty ticks" . 

The LM electronics can accumulate two-dimensional 
(2D) distributions of the multiplicity of in-time hits for 
the north and south LM detectors. An example distribu- 
tion, after background subtraction, accumulated at a lu- 
minosity of 63 fih~^s~^ is shown in Fig. [2] Three distinct 
components can be identified: (i) empty crossings with 
no LM hits, (ii) single-sided interactions where only one 
side has hits, and (iii) double-sided interactions where 
both sides have hits. 

In this study extensive use is made of these 2D mul- 
tiplicity distributions since they increase the number of 
events available for study by three orders of magnitude 
compared to an earlier analysis (Run Ila) and provide 
the ability to measure the multiplicity distributions for a 
single bunch crossing, instead of averaging over 36 bunch 
crossings. Consequently, rigorous background subtrac- 
tion techniques can be applied. In addition, data are now 
acquired over a short period of time (^8 mins total) to 
minimize the change in luminosity while the sample is ac- 
quired. For contrast, the Run Ila analysis includes ^1% 
statistical errors due to the measurement being based in 
low statistics data samples (on the order of ^ 10000 beam 
crossings). 
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Histogram data samples were acquired over a period 
ranging from August 2008 to January 2009 for a variety 
of luminosities. In total, 35 such datasets were used for 
this study. 



II. THE DO LUMINOSITY MEASUREMENT 

The DO luminosity measurement is performed by 
counting the rate of north-south coincidences in the LM 
detectors using 



L = 



1 dN 
a LAI dt 



(1) 



where ulm is the effective inelastic cross section seen by 
the LM, referred to as the "luminosity constant". The 
effective inelastic cross section is derived from the total 
inelastic cross section, UineU and adjusted for the LM 
system geometric acceptance and the efficiency for reg- 
istering inelastic events. The inelastic cross section has 
been measured at the Tevatron by the E710, E811, and 
CDF experiments Q. These experiments measure for- 
ward elastic scattering rates and use the optical theorem 
to determine the elastic, inelastic, and total pp cross sec- 
tion. A common averaging procedure [1] for the E811 and 
CDF measurements has been adopted by the CDF and 
DO experiments, which yields an inelastic cross section of 
ainei = 60.7 ± 2.4 mb at Vs = 1.96 TeV. 

The inelastic cross section can be subdivided into 
non-diffractive {und), single diffractive {(Jsd), and dou- 
ble diffractive (cdd) components. Single diffractive colli- 
sions are characterized by having the proton (antiproton) 
diffractively disassociate into hadrons while the antipro- 
ton (proton) remains intact. As in elastic collisions, the 
momentum transfer is typically small, so that the intact 
antiproton (proton) exits the detector through the beam 
pipe. Double diffractive collisions are similar to the sin- 
gle diffractive collisions except that both the proton and 
the antiproton undergo diffractive disassociation. Like 
single diffractive collisions, the particles produced tend 
to travel along the beam direction, thus producing large 
pseudorapidity gaps in the central region. Non-diffractive 
collisions represent the rest of the inelastic cross section 
and populate the full pseudorapidity region. Thus, the 
inelastic cross section can be expressed as 



fnd + f^sd + O'dd- 



(2) 



The effective inelastic cross section, aLM, can be writ- 
ten as 



CTLM — CFinel [fnd And + (1 ^ fnd) fsdAsd 

+ (1 - fnd) (1 - fsd) Am] 



(3) 



where fnd is the fraction of the inelastic cross section 
attributed to the non-diffractive process and fsd is the 



fraction of the diffractive cross section attributed to the 
single diffractive process, given by 



fnd — 
fsd 



f^nd 

7 

^inel 
(^sd 



<^dd 



(4) 
(5) 



The acceptances And, Asd, and Add are the non- 
diffractive, single diffractive, and double diffractive ac- 
ceptances, respectively, for producing at least one hit in 
both the north and south LM detectors. 

A single-sided effective cross section cjn (^s) can be 
defined for producing hits in only the north (south) LM 
detector 



llfndAnd + (1 fnd) fsdAg^ 

+ il-fnd){l~fsd)A^,] 



— <7inel[fndAnd + (1 ~ fnd) fsdAsd 

+ {1 - fnd) {1 - fsd) A^dl 



(6) 



(7) 



where A^j^\ Af^^\ A^j^^^ are the non-diffractive, sin- 
gle diffractive, and double diffractive acceptances, respec- 
tively, for producing at least one hit in the north (south) 
LM detector and no hits in the south (north) LM de- 
tector. Earlier analyses of the luminosity constant @ 
treated the north and south single-sided effective cross 
sections as being equal. In this analysis, small differences 
are found in both data and Monte Carlo (MC) simula- 
tions. These differences are attributed to asymmetries 
of the DO detector (e.g., the north endcap calorimeter is 
^ 4 cm closer to the pp interaction point than the south 
endcap calorimeter). Consequently, a separate calcula- 
tion of the north and south single-sided acceptances and 
cross sections is performed. 

Poisson statistics is used to determine the probability 
that a beam crossing will not have a north-south coinci- 
dence. This requires that there are no double-sided inter- 
actions characterized by (Tlm. In addition, there should 
be no pile-up of single-sided interactions that give hits on 
both the north and south sides stemming from multiple 
pp interactions in one beam crossing. Thus to contribute 
to the single-sided effective cross section, if there is a 
single-sided interaction on one side, there must be no 
single-sided interaction on the opposite side. 

The probability that a beam crossing has no north- 
south coincidences, P(0), is given by 



P(0) = 



(8) 



where L is the luminosity and / is the beam crossing 
frequency. The first factor is the probability for having 
no pp interactions giving a north-south coincidence. The 
term in the parenthesis corrects for two single-sided in- 
teractions in the same beam crossing mimicking a north- 
south coincidence. 
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The rate of live crossings with in-time hits in both the 
north and south LM detectors, Rlmi and the live cross- 
ing rate, Ruve, are measured in data. The probability 
for an empty beam crossing is derived from these rates 
to be 



P(0) = 1- 



Rl 



M 



Rlv 



(9) 



Given P(0), Eq. |5]is solved for the luminosity L, mak- 
ing use of the effective cross section ulm and the single- 
sided cross sections ctat and as- This calculation is per- 
formed for each of the 36 beam bunches independently 
since the luminosity, and thus P(0), is different for each 
bunch. 

This method is referred to as the "counting zeroes" 
technique. More details about this technique and its 
reach at high luminosity are presented in Appendix [Xj 



III. LUMINOSITY MONITOR BACKGROUNDS 

The LM is sensitive to two types of backgrounds: (i) 
out-of-time, and (ii) beam halo backgrounds. The out-of- 
time background is characterized by hits with an approx- 
imately uniform arrival time distribution for the ^40 ns 
measurement window before the beam crossing, with a 
significant variation in background rate over the 12 beam 
bunches in a bunch train. These hits randomly occur 
within the LM timing window, giving rise to in-time 
background hits. Studies of the out-of-time background 
have shown (see Appendix |B]) that the rate of these back- 
ground hits is proportional to the DO luminosity, indicat- 
ing that they are due to secondary particles from beam- 
beam interactions in previous beam crossings. This back- 
ground is found to have an effective single-sided cross sec- 
tion of 0.9 ±0.1 mb for both the north and south sides. It 
is assumed to originate mainly from low energy neutrons 
that interact in the LM scintillator, but it is possible 
that other sources, such as short-lived activation prod- 
ucts, contribute. 

Beam halo backgrounds occur when a proton or an- 
tiproton leaves the beam pipe upstream of the interac- 
tion point and produces secondary particles that are de- 
tected by the LM. Beam halo typically produces a shower 
in the upstream calorimeter that hits the upstream LM, 
and continues to the downstream LM. Typically, the up- 
stream LM counters will have out-of-time hits that arrive 
~9 ns earlier than particles from beam-beam collisions, 
while the downstream counters will have in-time hits that 
arrive at approximately the same time as particles from 
beam-beam collisions. Beam crossings with six or more 
early hits are vetoed in the luminosity calculation causing 
~1% of the beam crossings to be removed. The luminos- 
ity calculation is not sensitive to the selection of the veto 
value. 

During normal operation, the beam halo backgrounds 
are a few percent of the out-of-time background (e.g. for 



TABLE I: Background multiplicity rates distribution mea- 
sured in an empty tick immediately before a beam crossing. 
The data were acquired at a luminosity of 272 /ib~^s~^. The 
rows indicate the north multiplicity rate and the columns the 
south multiplicity rate, where the rates shown are normalized 
to the total number of events (5756304 events). 
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a luminosity of 100 /ib~^s~^ with a background cross sec- 
tion of 0.9 mb there is a 90 kHz background rate, whereas 
the typical halo rates are in the range of a few kHz) and 
have a negligible effect on the luminosity measurement. 
The remainder of this section focuses on describing and 
removing the out-of-time background. 

Table U shows a typical high luminosity background 
rates distribution measured in an empty tick immediately 
before a beam crossing. One or more background hits is 
present in 21% of these empty ticks, and 1.3% of them 
have hits in both north and south LM detectors. 

The largest contribution of the background in the de- 
termination of ulm occurs when an empty crossing is 
converted into a single-sided crossing with one back- 
ground hit. Previous determinations of the ulm for Run 
Ha Q have used two techniques for removing the back- 
ground contribution to the hit multiplicity distribution. 
The first technique used sidebands in the arrival time 
distribution to estimate the out-of-time background in 
events with a single hit in the detector and make the ap- 
propriate subtraction. The second technique required at 
least two hits (instead of one) on the opposite side when 
determining the non-diffractive fraction. 

For this study a different technique is used, for which a 
detailed description can be found in Appendix [C] If sig- 
nal and background hits are uncorrelated, the probability 
Dij of having i north counters and j south counters with 
observed hits is given by 

l<i p<i 
m<j q<j 

^ij = X! X! ^hnBpqflpifniqjQ {I + p - i) Q (m + q - j) , 

1=0 p=0 
m=0 q=Q 

(10) 

where Sim is the probability for having / north counters 
and m south counters with signal hits, Bpq is the prob- 
ability for having p north counters and q south counters 
with background hits, and Q is the Heaviside step func- 
tion 



e{l+p-i) 



1 ioT I + p > i 
ioT I + p < i. 



(11) 



The combinatoric factor fipi (and similarly fmqj) repre- 
sents the probability for observing i north counters with 
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FIG. 3: Example of the effect of background subtraction. The 
south multiphcity distribution is plotted for beam crossings 
where there are no hits in the north LM detector. The "Raw 
Data" histogram shows the multiplicity distribution before 
background subtraction. 



hits given that there were I north counters with signal 
hits and p north counters with background hits. The 
factor fipi is given by 



flpi 



ii{N~iy.p\ {N-p)\ 



(i+p- iy. (i - p)! (i - ly. (N - iy.m ' 



(12) 



where A'^ = 24 is the number of counters on a side. 

The observed signal and background probability distri- 
bution can be obtained from the multiplicity histogram 
associated with a beam crossing, and the background 
probability distribution can be obtained from the mul- 
tiplicity histogram obtained from the empty tick imme- 
diately before the beam crossing. Thus, we can solve 
the above set of linear equations for the background-free 
signal probability distribution. 

The effect of this background unfolding procedure is 
best illustrated by looking at its effect on beam cross- 
ings where one side has no hits. Figure [3] shows a slice 
of a 2D multiplicity distribution taken at a luminosity of 
63 nh^^s^^. The slice plotted shows the south multiplic- 
ity when there are no hits in the north LM detector. The 
prominent peak for a single south hit is substantially re- 
duced, as is expected from previous studies of the arrival 
time distribution for this bin. 

The background unfolding is performed on the 2D mul- 
tiplicity distribution, and it has certain features that do 
not appear in a ID unfolding procedure. For example, 
in Fig. [3] the number of entries in the higher multiplicity 
bins is larger after the background unfolding procedure 
than for the raw data. In the ID unfolding procedure, en- 
tries can only migrate from higher multiplicity to lower 
multiplicity. In a 2D unfolding procedure, entries can 
also migrate between ID slices. In the example shown 
in Fig. 131 diffractive events with no north signal hits but 
> north background hits are not included in the raw 



data (since the slice with no observed north hits is ex- 
amined), but are included in the background subtracted 
measurement since these events truly have no north sig- 
nal hits. The net effect of these two migrations is that 
bins with more than 3 south hits have more entries after 
background subtraction than before. 



IV. LUMINOSITY CONSTANT 
DETERMINATION 



A. Luminosity Monitor Acceptance Calculation 

The LM acceptances are calculated by simulating 
non-diffr active, single diffractive, and double diffractive 
events in the DO detector. Events are generated using 
PYTHIA [§1, with the CTEQ6L1 \Td\ parametrization of 
the parton distribution functions (PDFs), and utilizing 
the "Tune A" parameter set [ll| that is optimized to re- 
produce CDF data. The longitudinal distribution of the 
pp interaction vertex is generated as a Gaussian with an 
RMS of ~ 25 cm. In the transverse directions the beam 
is generated as a Gaussian centered at the origin with a 
width of about 100 /im. These Monte Carlo events are 
then processed through the standard DO detector simu- 
lation based on GEANtS [l^, using a modified detector 
geometry with the material model adjusted to match the 
multiplicity distribution in the LM as observed in the 
data. 

The number of hits that satisfy the timing criteria and 
have charge above threshold is counted separately for 
the north and south LM detectors, with the MC charge 
threshold and timing resolution adjusted to reproduce 
the data. Events are then classified into one of the follow- 
ing geometrical categories: (i) events that have at least 
one hit in both the north and south LM detectors (NS), 
(ii) events that have at least one hit in the north detector 
and no hits in the south detector (N Only), (iii) events 
that have at least one hit in the south detector and no 
hits in the north detector (S Only) and (iv) events that 
have no hits in either detector (Empty). The fraction of 
events in each category determines the associated accep- 
tance shown in Table HH The "N Only" and "S Only" 
acceptances are not identical since the DO detector, and 
its MC description, is not north/south symmetric. 



B. Determination of the Non-Diffractive Fraction 

The non-diffractive fraction can be related to the "zero 
fraction" measured in the data by using Poisson statis- 
tics. "Zero fraction" is the fraction of beam crossings 
that have no hits on a given side (north or south) when 
there is at least one hit on the opposite side. 

Starting with the north zero fraction, the probability 
of having no two-sided interactions and no one-sided in- 
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TABLE II: LM acceptances and their statistical uncertainties, for each type of inelastic process and geometrical category. 



Category 


Non-Diffractive 
Acceptance 


Single Diffractive 
Acceptance 


Double Diffractive 
Acceptance 


NS 


0.9924 ± 0.0009 


0.326 ±0.005 


0.500 ± 0.005 


N Only 


0.0048 ± 0.0007 


0.224 ±0.004 


0.203 ± 0.004 


S Only 


0.0026 ± 0.0005 


0.225 ±0.004 


0.212 ± 0.004 


Empty 


0.0002 ± 0.0001 


0.225 ±0.004 


0.0857 ± 0.0028 



teractions hitting the north side is 

P {N = 0) = e^^"^"'^"'^^^^ ^ . 



(13) 



Having zero hits on the north side while having at least 
one hit on the south side has a probability 



P (TV = 0, 5 > 0) = e-('"^^'+'"")^// [l - e-'^'^^lf^ 

Thus the north zero fraction is given by: 

^ _P{N = Q,S>0) _P{N = Q,S>Q) 



(14) 



Jo 



P{S>0) 
-(crLM+crN)L/ f ^ 



1-P{S 

-<ysL/f\ 



0) 



(15) 



1 _ Q-(TLM+<ys)L/ f 



The south zero fraction is obtained similarly, by exchang- 
ing (Tat with as- 

The cross sections in the above equations depend on 
the total inelastic cross section ainei, the non diffractive 
fraction /„d, the single diffractive fraction f^d, and the 
LM acceptances (see Eqs.[3l[6l[7l). The single diffractive 
fraction is taken to be fsd — 0.57 ± 0.21 |13j . 

Having already evaluated the LM acceptances, the only 
quantity still needed is the luminosity for the data sam- 
ple. In a given sample, the luminosity can be determined 
from the measured LM coincidence probability 



P (iV > 0, 5 > 0) = 1 - F (0) , 



(16) 



where P{0) is given by Eq. [51 Equations [12] and [1^1 are 
used to solve for the two remaining unknowns: the non 
diffractive fraction f„d, and the luminosity L. 

In the histogram data sample considered, 35 multiplic- 
ity histograms are acquired for the first beam bunch in 
the first beam train. For each of these histograms, the 
north and south zero fractions are calculated and the 
non-diffractive fraction determined. Results can be seen 
in Fig. [31 where the luminosity of the first beam bunch of 
the first beam train is extrapolated to the total delivered 
luminosity assuming that each beam bunch has the same 
luminosity. The statistical error on the non-diffractive 
fraction measurement ranges from 0.0006 to 0.0019 and 
is negligible. 

The north and south non-diffractive fractions are sim- 
ilar but exhibit some systematic differences. The north 
non-diffractive fraction has a higher dispersion, and ex- 
hibits a small correlation with luminosity which is at- 
tributed to the effects of beam halo, where the rate for 
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FIG. 4: Non-diffractive fraction measurements as a function 
of luminosity. Separate fits are performed to the north and 
south zero fractions. 



proton halo is in most cases substantially larger than for 
antiproton halo. Proton halo creates out-of-time hits in 
the north LM detector and in-time hits in the south LM 
detector. If a proton beam halo event occurs during an 
otherwise empty crossing, it will be counted as a single- 
sided event, thus increasing the north zero fraction. Since 
the diffractive processes are much more likely to have 
single-sided events than the non-diffractive processes, an 
increased zero fraction corresponds to a decrease in the 
non-diffractive fraction. The fraction of empty beam 
crossings decreases rapidly with increasing luminosity, so 
that the effect of beam halo is largest at low luminosity, 
which is consistent with the observed trend. The ob- 
served north/south difference is considered as a system- 
atic uncertainty in the calculation of (Tlm- Averaging 
the north and the south non-diffractive fractions for the 
35 luminosity points considered results in (/„d) = 0.668 
with an RMS spread of 0.002. 



C. Luminosity Constant 

The calculated LM acceptances f Table lll| and the non- 
diffractive fraction (evaluated in section ITVB[) . together 
with atnei = 60.7 ± 2.4 mb and fsd = 0.57 ± 0.21 are 
used to determine the luminosity constant, aLM, based 
on Eq.[31 for each data sample. Figure[51shows the consis- 
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FIG. 5: Measured luminosity constant as a function of time. 



FIG. 7: Distribution of luminosity constant measurements. 
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FIG. 6: Measured luminosity constant as a function of the DO 
luminosity. 



tency of the measurement of (Tlm during the period of «6 
months when the muhiplicity histograms were collected. 
The luminosity constant is labeled as "North" when the 
north non-diffractive fraction has been used, while it is 
labeled as "South" when the south non-diffractive frac- 
tion has been used. Figure [S] shows the luminosity de- 
pendence of the CTLAf determination, which indicates a 
similar trend to the one observed for the non-diffractive 
fraction (see Fig. |4]). The statistical uncertainty on the 
clm measurements ranges from 0.02 to 0.07 mb. 

The distribution of the ctlm measurements is shown in 
Fig. [71 The average of the north and south luminosity 
constants gives 



o'LM ~ 48.3 mb. 



(17) 



The north and south single sided cross sections are de- 
termined to be 



(Tat 



as = 4.5 mb. 



(18) 



D. Luminosity Constant Uncertainty 

In evaluating the systematic uncertainties affecting the 
measurement of ctlm, the different sources are propa- 
gated through the analysis chain to establish the effect 
on <jlm, including re-calculating the LM acceptances and 
determining the non-diffractive fraction for each data 
sample. This allows us to take into account the correla- 
tion between the LM acceptances and the non-diffractive 
fraction. The different sources of systematic uncertain- 
ties are listed below: 

Inelastic Cross Section 

The CDF and DO experiments have adopted [3| an 
inelastic cross section of Uinei — 60.7 ± 2.4 mb at 
^/s — 1.96 TeV for their luminosity measurements. 
Propagating the 2.4 mb uncertainty in the inelastic 
cross section gives an uncertainty of 1.91 mb on 
the luminosity constant, (Tlm- 

Single DifFractive Fraction 

The single diffractive fraction is taken to be 
fsd — 0.57 ± 0.21 [l3|, corresponding to a large 
variation in the single diffractive cross section. 
This corresponds to an uncertainty of 0.43 mb on 

Time Variation / Radiation Damage 

Periodic adjustments to the PMT high volt- 
age are performed, and the LM scintillators are 
replaced during long shutdowns to minimize the 
impact of radiation damage. The high voltage 
changes typically lead to less than 0.5% change in 
the measured luminosity. An uncertainty of ±0.5% 
(0.24 mb) is assigned to glm due to time variation 
in the luminosity measurement. 
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Hit Multiplicity 

FIG. 8: Hit multiplicity when there are one or more hits on 
the opposite side. The north and south distributions have 
been combined to reduce MC statistical errors. The data are 
from a background subtracted multiplicity histogram taken 



Hit l\Aultiplicity 

FIG. 9: Hit multiplicity when there are no hits on the oppo- 
site side. The north and south distributions have been com- 
bined to reduce MC statistical errors. The data are from a 
background subtracted multiplicity histogram taken at a lu- 
at a luminosity of 63 ^b~^s~^. The MC samples represent minosity of 63 /xb^^s^^. The MC samples represent the tuned 



the tuned material model. 



material model. 



GEANT Energy CutofTs 

The simulated events are reprocessed through 
the DO GEANT3 simulation with lower energy 
cutoffs. The S ray production cutoff is lowered 
from 1 MeV to 10 keV, the neutral and charged 
hadron cutoffs are lowered from 1 McV to 100 kcV, 
and the muon cutoff is lowered from 10 MeV to 
100 keV. The 100 keV energy cutoff selected for 
this study is based on an estimate of the lowest 
energy that could cause an LM counter to detect a 
hit. A change of 0.24 mb in itlm is observed, and 
this is assigned as a systematic uncertainty. 

Monte Carlo Material Model 

In order to have the MC reproduce the ob- 
served hit multiplicity distribution, it is necessary 
to modify the material in front of the LM in the 
GEANtS model. MC events have been generated 
both with the nominal (untuned) material model 
and the tuned material model. The resulting 
change in aLM of 0.16 mb is taken as a systematic 
uncertainty. 

As a check that the tuned material model repro- 
duces the hit multiplicity histograms, Figs. [S] and [S] 
show the hit multiplicity for a histogram acquired 
at a luminosity of 63 fib~^s~^. The data points 
have the out-of-time background contribution re- 
moved using the background unfolding procedure 
described in Appendix [C] In Fig. |8] at least one 
opposite side hit is required, while in Fig. |9] the 
requirement is that there are no opposite side hits. 
The hit multiplicity distributions for these data and 



MC samples are observed to be in good agreement. 

Luminosity Monitor Acceptance 

The acceptance of the LM shows a small de- 
pendence (less than 0.1%) due to the longitudinal 
variation in the location of the pp interaction 
vertex. In addition, the variations in the trans- 
verse position of the interactions^ results in a 
0.17% variation in the LM acceptance. Adding 
these effects in quadrature with the Monte Carlo 
statistics yields an uncertainty of 0.11 mb on glm- 

Light Collection / Radiation Damage 

Radiation daniage can reduce the light collec- 
tion efficiency [Tl. To estimate this contribution, 
MC events with a piece-wise linear change in the 
light collection efficiency were generated. The 
light collection efficiency is reduced by a factor of 
two at the inner edges of the scintillator wedge, 
increasing linearly to no change at the center of 
the PMT, and then decreasing linearly to a factor 
of two reduction at the outer edge of the wedge. 
The charge threshold is also adjusted to emulate 
the effect of the HV changes that work to keep 



^ The beam position could vary as much as 0.3 mm around the 
nominal transverse location during normal data taking, and the 
nominal transverse location was 1.6 mm off center before October 
2007 and was moved to 0.2 mm off center for stores after that 
time. 
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the average charge constant. These changes in the 
MC efficiency calculation resulted in a change in 
of 0.09 mb, which is assigned as a systematic 
uncertainty. 

North — South Asymmetry 

The north and south (Jlm measurements dif- 
fer by an average of 0.18 mb. Since the final value 
of (Jlm is selected to be the average of the north 
and south measurements, a systematic uncertainty 
of 0.09 mb is assigned to aLM to account for the 
observed north - south asymmetry. 

Luminosity Dependence 

In the ensemble of 35 <Jlm measurements 
made at varying luminosities, an RMS spread of 
0.08 mb is observed, so this value is assigned as a 
systematic to account for the observed luminosity 
dependence. 

PDF Choice 

The standard DO PDF set is CTEQ6L1 Events 
were also generated using MRST2004NLO [H]. 
This change resulted in a difference of 0.06 mb in 
clm, which is assigned as an uncertainty. 

PYTHIA Tune 

Events have been generated with a modified 
value for the transition point between the low-px 
and high-pT models, as an alternative to modifying 
the material model. This is done by changing 
the PYTHIA 9] parameter PARP(82), which is 
known to have a significant affect on the average 
multiplicity, from the default value of 2.0 to 1.25. 
This change in the pythia parameter led to a 
change in aLM of 0.04 mb, which is assigned as a 
systematic uncertainty. 

Background Unfolding 

The (Tlm calculation is repeated using multi- 
plicity histograms acquired during the last bunch 
crossing of a bunch train. These bunches have 
«40% higher background than the first bunch of 
a bunch train. The difference of 0.03 mb between 
the two calculations is assigned as a system- 
atic uncertainty associated with the background 
unfolding. 

GEANT Hadronic Model 

For this study, the GCalor I16t hadronic model 
is replaced with the Geisha [l3| model due to 
its ability to better handle low energy particle 
interactions. This change in the GEANT hadronic 
model resulted in a 0.03 mb change in glm, which 
is assigned as a systematic uncertainty. 



TABLE III: Contributions to the luminosity constant {(jlm) 
uncertainty. 

Source Uncertainty (mb) 



Seasonal Timing Variation 

There are seasonal drifts in the DO clock stem- 
ming from temperature variations that result in 
expansion or contraction of the long cable used 
to send signal for collisions from the accelerator 
control room. To account for the effect of seasonal 
timing observed, the timing window of ±6.4 ns 
that defines an in-time hit is shifted by ±1 ns. The 
resulting change of 0.02 mb in glm is assigned as 
a systematic uncertainty. 

Charge Threshold 

The charge threshold in the MC simulation is 
shifted by ±2 pC around the nominal value of 8 pC 
to account for uncertainties in the modeling of the 
charge threshold. The resulting change of 0.01 mb 
in ulm is assigned as a systematic uncertainty. 

The luminosity constant {(Jlm) uncertainties are 
summarized in Table IIIII Adding the uncertainties in 
quadrature yields an uncertainty in ulm of ±2.0 mb, 
where 1.9 mb is associated with the uncertainty in the 
inelastic cross section, 0.4 mb is associated with the 
uncertainty in the single diffractive fraction, and 0.4 mb 
is associated with the remaining sources of uncertainty. 



V. INTEGRATED LUMINOSITY 

The integrated luminosity £, defined as 

L = I L-dt, (19) 

where T is the data taking period, is the relevant quantity 
used in the measurements of cross sections and in setting 
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FIG. 10: The recorded luminosity profile in bins of 5 jJo ^s 



upper limits on the production of new particles. In ad- 
dition to the uncertainty on the luminosity constant, the 
uncertainty on the determination of the integrated lumi- 
nosity takes into account additional contributions that 
cover possible variations with time and with luminosity 
of the luminosity constant. These additional sources of 
uncertainty are discussed below. 

The "delivered" luminosity is the integrated luminosity 
delivered by the Tevatron. The "recorded" luminosity is 
the integrated luminosity associated with a specific trig- 
ger and takes into account the deadtime and losses in the 
data acquisition system. Level I triggers are grouped to- 
gether so that they have common deadtime, i.e., common 
sources of enable, disable, and readout The recorded 
integrated luminosity referred to in this section corre- 
sponds to the luminosity exposure of the experiment's jet 
trigger with the highest transverse energy, = E ■ sin0, 
which requires at least one jet with Et > 125 GeV. 

The total recorded integrated luminosity for Run lib is 
assessed to be 9.2 fb~^. Of that, ^0.2 fb^^ were recorded 
with luminosity above 300 /ib^^s"^, and ^-^1.7 fb~^ were 
recorded with luminosity between 200 and 300 /ib^^s"^. 
Figure [10] shows the recorded luminosity profile, in bins 
of width 5 /ib^^s^^, for the full Run lib dataset. 



Long-Term Stability of the Luminosity 
Measurement 



The yield of single muons in the forward muon sys- 
tem [13] can be used as an independent check of the sta- 
bility of the luminosity measurements. The single muon 
yield Y is monitored regularly using special data samples 
and is obtained by 



Y 



N, 



(20) 



where is the number of the forward muons, and C is 
the integrated luminosity of the respective data sample. 
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FIG. 11: Normalized single muon yields (a) 
time, and (b) its projection on the y-axis. 



as a function of 



The normalized yield is obtained by dividing the value of 
the muon yield for each data sample by the mean value 
of all the data samples collected 



Y — 



Y 



mean value 



(21) 



Figure 11(a) shows the normalized yield measurements 



for Run lib. If the statistical error on TV^ is less than 1% 
the uncertainty on the yield, cry, is set to 1%; otherwise 
cry is set to the statistical error of N^. Figure 11(b) 
shows the distribution of the normalized muon yields 
super-imposed with a Gaussian function that indicates 
a typical variation of ~0.8%. 

In addition, muon yields are measured as a function 
of luminosity. Figure [12] shows an example of the dis- 
tribution of the normalized muon yields collected during 
two Tevatron stores. The super-imposed Gaussian on 
Fig. [12(b)] indicates a typical variation of ~0.4%. After 
accounting for the statistical errors on A'^ and the un- 
certainty due to the time variation / radiation damage 
already included in the luminosity constant, a systematic 
uncertainty of 0.6% is assigned to the integrated luminos- 
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FIG. 12: Normalized single muon yields (a) as 
luminosity, and (b) its projection on the y-axis. 



function of 



ity measurement based on the observed variation in the 
muon yield. 



B. Integrated Luminosity Measurement 
Uncertainty 

The uncertainty of the integrated luminosity includes 
a dominant contribution of 4.2% stemming from the un- 
certainty on the luminosity constant. Additionally, it 
includes an uncorrelated contribution of 0.6% in the inte- 
grated luminosity measurement described in the previous 
section. Other potential source of systematic uncertainty 
includes events that are lost in the DO data acquisition 
and reconstruction system and are not properly taken 
into account. The fraction of these events has been stud- 
ied and found to be negligible (<< 0.1%). 

Adding these uncertainties in quadrature yields an un- 



certainty in the integrated luminosity of 4.3%, where 
4.0% originates from the uncertainty in the inelastic cross 
section and is correlated with both the CDF and DO Run 
Ila integrated luminosity measurements, 0.9% originates 
from the uncertainty in the single diffractive fraction and 
uncorrelated with the CDF integrated luminosity mea- 
surement, but correlated with the DO Run Ila integrated 
luminosity measurement, and 1.0% originates from the 
remaining sources of uncertainty that are uncorrelated 
with both the CDF and DO Run Ila integrated luminos- 
ity measurements. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have measured the effective inelas- 
tic cross section, aLM, as seen by the DO luminosity 
monitor, and assessed the recorded integrated luminos- 
ity for data collected with the DO detector at the Fer- 
milab Tevatron Collider for the period called Run lib 
(June 2006 to September 2011). A luminosity constant 
of aLM — 48.3 ± 2.0 mb is obtained. The Run lib lumi- 
nosity constant is 0.6% larger than the Run Ila luminos- 
ity constant (48.0 ±2.9 mb) and its uncertainty has been 
reduced from 6.1% to 4.2%. 

The recorded integrated luminosity for the highest Et 
jet trigger is £ = 9.2 ± 0.4 fb"! during Run lib. The 
total relative uncertainty of the Run lib DO recorded in- 
tegrated luminosity is determined to be 4.3%, where 4.0% 
is associated with the inelastic cross section, 0.9% is as- 
sociated with the single diffractive fraction, and 1.0% is 
associated with the DO sources of uncertainty. 
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Appendix A: Empty Crossing Method 

The DO luminosity is derived from the number of beam 
crossings with north - south in-time coincidences in the 
DO LM that occur during a measurement period. The 
luminosity reported to the Fermilab accelerator division 
for monitoring purposes is based on a nominal 15 s mea- 
surement period, while the luminosity used for physics 
analyses employs a nominal 60 s measurement period. 
The actual measurement period is occasionally shorter 
to ensure synchronization between the luminosity mea- 
surement and the state of the data acquisition system. 
Statistical fluctuations in the number of luminosity coin- 
cidences lead to statistical and systematic errors in indi- 
vidual luminosity measurements. 

The luminosity calculation is performed separately for 
each of the 36 beam bunches in the Tevatron, resulting 
in each bunch having its own measured luminosity. For 
a bunch with true luminosity L, the average number of 
proton - antiproton interactions that produce north - 
south coincidences is proportional to the luminosity 



(Al) 



where ctlm is the effective inelastic cross section seen by 
the LM and / is the beam crossing frequency. Similarly, 



the average number of interactions in a beam crossing 
that produce in-time hits in only the north (south) LM 
array is given by 



Ns{L) = ^L, 



(A2) 



where crjv and as are the single-sided efii'ective cross sec- 
tions. 

The fraction of beam crossings that do not produce a 
north-south coincidence of in-time LM hits are classified 
as empty and, using Poisson statistics, is given by 

Fo(L) = e-^^"(^) 

(A3) 

where the first factor is the probability for having no pp 
interactions giving a north-south coincidence. The term 
in the parenthesis gives the probability for not having 
multiple single-sided interactions that result in a north- 
south coincidence. 

The average number of empty crossings during a mea- 
surement period is then 



No{L) = NuveXFoiL), 



(A4) 



where A^Live is the number of live beam crossings during 
the measurement period. In the results shown below, 
we ignore the small fraction of beam crossings that are 
rejected due to halo veto ( -^1%) and we consider A^Live = 
/ • T, where T is the measurement period. 



1. Behavior of the DO Luminosity Measurement at 
High Luminosity 

The number of empty crossings observed in difi^erent 
measurement periods can be described with the bino- 
mial distribution. In the high luminosity limit where the 
number of empty crossings is small, the distribution of 
the observed number of empty crossings is well approxi- 
mated by a Poisson distribution 



Pino) 



,Wo(L) ^ 



-No{L) 



(A5) 



where P(no) is the probability of observing no empty 
crossings. 

For each measurement period, the measured luminos- 
ity Lm corresponding to the observed number of empty 
beam crossings, no, is calculated by numerically solving 
the equation 



no = ATii^ee-'^^"^'"/^ 



(A6) 
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In the case where no empty crossings are observed, the 
solution of the above equation yields an infinite measured 
luminosity. When that occurs, the luminosity is set to the 
value that would be found if there had been one empty 
crossing observed. The impact of this approximation on 
the measurement of the luminosity in DO is discussed 
below. 

The average measured luminosity is given by 



500 



P(0)L„(1) 



fT 

"■0 = 1 



(A7) 



where P{no) is the probability of observing uq empty 
crossings, and the first term accounts for the special han- 
dling where no empty crossings are observed. 

The total luminosity is obtained by summing the lu- 
minosity from the 36 beam bunches. If all 36 bunches 
had the same luminosity, the total luminosity would be 
36 times the bunch luminosity and the RMS spread of 
the total luminosity would be a factor of six times the 
RMS for a single bunch since each bunch measurement is 
statistically independent. In practice there are typically 
a few percent variations among the bunch luminosities. 
While these bunch-to-bunch variations are accounted for 
in the DO luminosity measurement, the results following 
illustrate the behavior of the measured total luminos- 
ity under the assumption that all bunches have the same 
bunch luminosity. Including typical bunch-to-bunch vari- 
ations will not significantly affect these results. 

The relation between the average measured luminos- 
ity and the true luminosity for measurements periods of 
15 s and 60 s is shown in Fig. [131 Due to the non-linear 
behavior in the empty crossing probability at high lumi- 
nosity, the average measured luminosity systematically 
exceeds the true luminosity before entering the saturation 
region, where the luminosity asymptotically approaches 
the value L„i{l) that is assigned to bunches with less 
than two empty beam crossings. 

Table ITVl illustrates how the non-linear behavior of the 
empty crossing probability and the special treatment of 
the case where there are no empty crossings in the mea- 
surement period lead to non-linear behavior for the av- 
erage measured luminosity. In this example, the true 
luminosity is 420 /ib~^s~^, the measurement period is 
15 s, and the average number of empty crossings in the 
measurement period is 3.05. The exponential decrease 
in the measured luminosity as the number of observed 
empty crossings increases results in an average measured 
luminosity of 424.8 /ib~^s~^, which is 1.1% higher than 
the true luminosity. 

Figure [14] shows the mean deviation between the av- 
erage measured luminosity and the true luminosity as a 
function of the true luminosity. The largest positive devi- 
ation occurs when the average number of empty crossings 
is ~ 3. This occurs at a luminosity of 420 ^b~^s~^ for 
a 15 s measurement period, where the average number 
of interactions with a north-south coincidence is ^ 12 
and the fraction of empty crossings is ~ 4 • 10~®. For a 
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FIG. 13: Average measured luminosity versus the true lumi- 
nosity for 15 s and 60 s measurement periods. 



TABLE IV: Non-linearity in the luminosity measurement 
showing the probability, measured luminosity, and contribu- 
tion to the average luminosity as a function of the observed 
number of empty crossings. In this example, the true lumi- 
nosity is 420 /ib~^s"^ and the measurement period is 15 s. 



no 


P{no) 


Lm (^tb^^s"^) 


P{no)-Lm {^ih-^s-^) 





0.047 


456.9 


21.7 


1 


0.145 


456.9 


66.1 


2 


0.220 


434.0 


95.7 


3 


0.224 


420.5 


94.2 


4 


0.171 


411.0 


70.1 


5 


0.104 


403.6 


42.0 


6 


0.053 


397.5 


21.0 


> 6 


0.036 




14.0 


Sum 


1 




424.8 



60 s measurement period, an average of 3 empty cross- 
ings occurs at a luminosity of 470 ^h^^s~^, where the 
average number of interactions with a north-south coin- 
cidence is ~ 13 and the fraction of empty crossings is 
~ 10~^. Given that only a small fraction of data was 
collected at luminosities in excess of 300 ^b~^s~^ and 
that the measurement period used for the determination 
of the integrated luminosity is of 60 s, we estimate that 
the impact of the non-linear behavior in the empty cross- 
ing probabilities has a negligible impact on the precision 
of the integrated luminosity determination in DO. 

The RMS width of the measured luminosity for a single 
bunch is given by 



fT 



"0 = 1 

(A8) 

where the first term accounts for the special handling 
where there are no empty crossings observed. The sta- 
tistical uncertainty in the total luminosity for the 36 
bunches is shown in Fig. 1151 The statistical uncertainty 
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FIG. 14: Mean deviation between the average measured lumi- 
nosity and the true luminosity for 15 s and 60 s measurement 
periods. 
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FIG. 15: RMS width of the measured luminosity for 15 s and 
60 s measurement periods. 



in the luminosity measurement is less than 0.1% for lumi- 
nosities in the range 1.5 ~ 250 fj,h~^s~^ for the 15 s mea- 
surement period and 0.4 ~ 310 fjh~^s^^ for the 60 s mea- 
surement period. The RMS width reaches a maximum 
of 0.7% when the average number of empty crossings 
is ~ 3. Further increases in the true luminosity push the 
luminosity measurement into the saturation region where 
an increasing fraction of luminosity measurements report 
the maximum possible measured luminosity Lm(l), lead- 
ing to a decrease in the RMS width. 

At high luminosity, the statistical error can be approx- 
imated by 



L 



1 



1 



/NbNns Ro 



(A9) 



where i?o is the rate of empty crossings, Ng is the num- 
ber of beam bunches, and the small contribution to the 



empty crossing rate from multiple single-sided interac- 
tions is ignored. Thus, the large statistical error on the 
empty crossing rate when there are an average of 3 empty 
crossings {l/VS or 58%) is reduced by a factor of six due 
to the 36 independent bunch measurements and an ad- 
ditional factor of ^ 12 at 420 fih^^s^^ due to the expo- 
nentially falling empty crossing rate to yield a precision 
in the total luminosity of better than 1%. 

A key requirement for the empty crossing method to 
work at high luminosity is that beam crossings with in- 
time north-south coincidences are not misclassified as 
empty crossings. With empty crossings rates in the part- 
per-million range at peak luminosities, this misclassifi- 
cation probability must be well under 10^^ for typical 
beam crossings. Such beam crossings are easily identi- 
fied in the DO LM since a typical beam crossing yields 
in-time hits in most or all of the 24 luminosity counters, 
whereas the requirement is > 1 in-time hit. Since early 
hits in a large number of counters can mask the pres- 
ence of in-time hits, a halo veto has been implemented 
to exclude beam crossings with more than six early hits 
from the luminosity calculation. Timing distributions are 
monitored and the LM TDCs are recalibrated as needed 
to ensure that the timing distributions are well centered 
within the timing window. Tests of the digital counting 
electronics have proven this system to be robust with no 
evidence of misclassification. 



Appendix B: Background Studies 
1. Luminosity Dependence of Background Rate 

The luminosity dependence of the rate for back- 
ground hits in bunch crossings that do not have proton- 
antiproton collisions is studied, where the multiplicity 
of background hits is neglected and only the fraction of 
beam crossings with > hits in the LM is measured. The 
north and south background rates, and the coincidence 
rate where background hits occur in both north and south 
counters are separately measured. Figure 16(a) shows the 
background rates as a function of luminosity, measured in 
an empty tick 132 ns before the first bunch crossing of a 
bunch train. Figure [l 6 (b) | shows the same rates measured 
in an empty tick 396 ns after the last bunch crossing of 
a bunch train. In both cases, the north and the south 
rates scale approximately linearly with luminosity, while 
the north - south coincidence rate is much lower and has 
a non-linear dependence on luminosity. 



2. Effective Background Cross Section 

Since the north and south background rates scale lin- 
early with the luminosity, the background rate can be 
treated as an effective background cross section. Since 
there is no actual luminosity in the ticks where the back- 
grounds are measured, 1/36 of the DO luminosity (i.e.. 
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the average luminosity attributable to one of the 36 beam 
bunches) is taken to calculate the background cross sec- 
tion. 

To account for pileup in the background, Poisson 
statistics is used to relate the background cross section 
and the background rate Rbg in a given detector using 

Rbg 



P(0) = e"""^^^/^ = 1 - 



/ 



f , /, Rbg 



(Bl) 
(B2) 



where aBG is the background cross section, L is the lu- 
minosity, / is the beam crossing frequency, and Rbg is 
the background rate. 

Figure [17] shows that the background cross section is 
largely independent of the luminosity. The background 
cross section after a bunch train is about 40% higher 
than immediately before the bunch train, indicating that 
the background cross section increases during the bunch 
train. 

For a given bunch crossing, the background cross sec- 
tion is observed to decrease over time, as shown in 
Fig. [iHl and attributed to radiation damage to the scin- 
tillator. The cross section increases significantly in data 
recorded in late January 2009 immediately following an 
increase in the high voltage applied to the LM PMTs to 
compensate for the radiation damage. There is also an 
increase in the background rate in early October 2008, 
following a long downtime of the accelerator, interpreted 
as evidence of some annealing of the radiation damage 
in the scintillator following an extended shutdown. Av- 
eraging the 140 background cross section measurements 
from 35 data samples shown in Fig. [171 the average back- 
ground cross section is estimated to be: asG = 0-85 mb. 



TABLE V: Background multiplicity distribution measured in 
an empty tick immediately before a beam crossing. The data 
were acquired at a luminosity of 272 ^b~^s~^. The rows in- 
dicate the north multiplicity and the columns the south mul- 
tiplicity. 



N/S 
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> 4 





4522672 


466907 
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3203 
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501848 


52260 
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394 
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84467 


9072 


1524 


208 


63 
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11706 


1349 


263 


33 


17 


> 4 


3460 


399 


78 


23 


10 



i counters hit in the north luminosity monitor and Ps{j) 
is the probability of having j counters hit in the south 
luminosity monitor. 

If the backgrounds in the north and the south detec- 
tors are uncorrelated, the number of entries in the 2D 
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3. Background Characteristics 

The coincidence background where both north and 
south luminosity monitors recorded hits is studied, since 
such coincidences are counted in the luminosity measure- 
ment. Table fVl shows a high luminosity background mul- 
tiplicity distribution for an empty tick immediately be- 
fore a bunch train. One or more background hits are 
present in 21% of the beam crossings, and 1.3% of the 
beam crossings have hits in both north and south lumi- 
nosity monitors. The multiplicity of background hits is 
low. Approximately 79% of beam crossings with back- 
ground hits have only one background hit among the 48 
luminosity counters. 

The question whether the north and south background 
hits are statistically independent is studied by projecting 
the 2D distributions to obtain ID probability distribu- 
tions for observing a particular number of hits in the 
north/south monitor. Table IVll shows the ID probability 
distributions derived from the 2D multiplicity distribu- 
tion in Table fVl where PAr(i) is the probability of having 
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FIG. 16: Background rate, measured in an empty tick (a) 
immediately before the first bunch crossing and (b) 396 ns 
after the last bunch crossing of a bunch train. One or more 
background hits in the north/south luminosity monitors is 
required (left axis). Also shown is the rate of beam crossings 
in which both north and south luminosity monitors have one 
or more background hits (right axis). 
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TABLE VI: North and south ID probabihty distributions de- 
rived from the 2D multiphcity distribution in Table [V] 



FIG. 17: Effective background cross section as a function of 
luminosity, measured in an empty tick 132 ns before (tick 6) 
and 396 ns after (tick 43) a bunch train. 
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FIG. 18: Time dependence of the background cross sections. 
The increase in the cross section in late January 2009 corre- 
sponds to the raising of the PMT high voltages to compensate 
for radiation damage. The PMT high voltages were adjusted 
a few days before the first data point in the plot. Tick 6 refers 
to an empty tick immediately before a bunch train, and tick 
43 refers to an empty tick after a bunch train. 



multiplicity distribution, N{i,j), will be given by 

Nit,j)^PN{i)Ps{j)No, 



(B3) 



where A^o is the total number of entries in the 2D mul- 
tiplicity distribution. Table IVIII shows the predicted 2D 
multiplicity distribution using this equation and the ID 
probability distributions in Table IVII Reasonably good 
agreement is observed between the measured 2D multi- 
plicity distributions shown in Table |V] and the predicted 
distribution in Table IVIII Nevertheless, there is a small 
but persistent underestimate in the number of predicted 
coincidences between north and south detectors. 
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FIG. 19: Difference between the measured and predicted rates 
before (Tick 6) and after (Tick 43) bunch train crossings with 
hits in both north and south luminosity monitors. 



Figure [19] shows the difference between the predicted 
and measured rate for observing north - south back- 
ground coincidences for the data sample considered. The 
predicted rate is consistently underestimated by a small 
amount, and the magnitude of the discrepancy grows 
with luminosity. 

There are no additional known sources of background 
that can give north - south coincidences. In predicting 
the random coincidence rate, the assumption was made 
that there are no correlations in the probabilities that the 
north and south detectors have background hits. Since 
the rate of background hits scales with luminosity, the 
background must be associated with beam ~ beam in- 
teractions in previous beam crossings. The number of 
interactions in a given beam crossing will fluctuate ac- 
cording to Poisson statistics. Beam crossings with up- 
ward fluctuations in the number of interactions will have 
a higher probability of producing background hits in both 
the north and south luminosity monitors, while crossings 
with a downward fluctuation will have a lower probability 
of producing background hits. 

Let P represent the probability of producing a back- 
ground hit in either the north or south luminosity moni- 
tor for a given beam crossing. Instead of assuming a fixed 
value for P, it is assumed that P fluctuates depending 
on how many beam - beam interactions took place in 
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TABLE VII: Predicted 2D probability distribution obtained 
from the ID probability distributions in Table IVII under the 
assumption that the north and south ID distributions are 
statistically independent. 
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recent beam crossings. The probability of having a north 
- south coincidence, Pcoin, is given by 



(Pn. 



= in 



(B4) 



where ap is the RMS spread of the P distribution. If 
there is an average of -/Ve// beam - beam interactions 
producing background hits, we obtain: 



and 



(P) 



= (P) 



(B5) 



(B6) 



Assuming that each of the 36 beam crossings has the 
same luminosity, L — Ltot/iQ, and that the background 
hits effectively arise from the previous m beam crossings, 
we obtain 



(B7) 



The parameter m is estimated by fitting the data. Us- 
ing empty ticks prior to a bunch train crossing gives an es- 
timate of m « 7 and after the bunch train crossing m w 9. 
The difference between the predicted and measured north 
- south coincidence rates, both with and without correct- 
ing for fluctuations in the background probability, are 
shown in Fig. [20l The background model gives good 
agreement between the predicted and measured coinci- 
dence rate provided that the background probability is 
allowed to fluctuate. 



4. Background during Beam Crossings 

The aim is to estimate the background cross section 
for beam crossings, but the method discussed up to now 
is not applicable to actual beam crossings. A model is 
developed that allows the estimation of the background 
cross section during beam crossings. Using this model, it 
is found that the average background cross section for the 
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FIG. 20: Difference between the measured and predicted rates 
for bunch train crossings with hits in both north and south 
luminosity monitors (a) before and (b) after the bunch train. 
The prediction of Eq. IB3I assumes a fixed background proba- 
bility, whereas the corrected prediction allows the background 
probability to fluctuate. 



36 beam crossings is slightly higher than what is obtained 
by averaging the background cross sections before and 
after the beam crossing. 

The probability of producing a background hit in the 
north or south luminosity monitors is estimated from the 
coincidence probability 



P = 




(B8) 



The corresponding background cross section is calcu- 
lated using Poisson statistics to account for pileup 



-^ln(l-P), 



(B9) 



Figure m] shows the background cross section for data 
recorded over the span of one minute on November 12, 
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FIG. 21: Background cross section during bunch crossings. 



FIG. 23: Predicted background cross section. 



2008 at a luminosity of 68 /xb^^s^^. The background 
cross section rises during the course of bunch crossings 
and then decays in the region without bunch crossings. 
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FIG. 22: Comparison of the predicted background cross sec- 
tion and the measured background cross sections. 

A model is constructed to describe the observed back- 
ground cross section. The amplitude of the background 
produced by a given beam bunch is taken to be propor- 
tional to the luminosity for that bunch. It is assumed 
that the contribution to the background cross section 
falls exponentially in time following the bunch. There 
appears to be both a short and long time components in 
the background, thus a double exponential is used to fit 
the background. For a bunch occurring at t = with 
luminosity L, the background model predicts that the 
contribution of the background to future crossings will 
be 



o-Bcit) = L ( Aie 



A2e 



-t/T2 



(BIO) 



where Ai and A2 are the amplitudes of the two back- 



ground components and ti and T2 are the associated time 
constants. The background cross sections are fitted us- 
ing this model. As shown in Fig. 1221 this four parameter, 
two-component model provides a good description of the 
data. In this particular example, the gap region shows an 
anomalous up-tick in the background cross section that 
is likely due to residual beam in these ticks. In fitting the 
background model, the three ticks with more prominent 
anomalies have been excluded from the fit as well as the 
two preceding ticks. 

Since this model describes adequately the tick depen- 
dence of the background cross section, it is also used to 
estimate the background cross section during the beam 
crossings (see Fig. [23|) . The background cross section 
varies over time, thus the final background cross sec- 
tion estimate is based on the average over the 35 data 
samples used for this study. Previously (see Sec. IB 21) . 
a background cross section of asG = 0.85 mb was es- 
timated by averaging the individual measurements for 
north and south monitors, before and after beam cross- 
ings. Correcting this value according to the model de- 
scribed above, a background cross section of asG — 
0.9 ± 0.1 mb is calculated, where the uncertainty quoted 
is the RMS spread. 



Appendix C: Background Subtraction 

The analytical method used to disentangle the back- 
ground contribution from the measured multiplicity dis- 
tributions is described below. 



1. Signal and Background Convolution 

Given i counters with hits in one of the LM arrays, lim- 
its can be placed on the number of signal and background 
hits even though the same configuration of hits can be 
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obtained from different combinations of signal and back- 
ground hits. If I counters have signal hits and p counters 
have background hits, I and p will be constrained hy I < i, 
p <i, and I +p> i. The last constraint is not an equal- 
ity because there can be counters with both signal and 
background hits in them as multiple particles hitting a 
given counter are counted as a single hit. 

The probability of having i counters with hits can be 
constructed if the multiplicity distributions for signal and 
background are known. If Si and Bp are the signal and 
background probabilities, then the probability dup for 
having i observed counters with hits is given by 



Hip 



SlBpflpi , 



(CI) 



where fipi is a combinatoric factor that gives the proba- 
bility that I counters with signal hits and p counters with 
background hits will yield i counters with observed hits. 
The combinatoric factor fipi is derived in Section IC 21 

The total probability Di for observing i hits is obtained 
by summing over the possible values of I and p: 

D^^J2Y. S^Bpflp^Q {I + P - i) , (C2) 

1=0 P=0 

where the constraints are explicitly represented through 
the summation limits and the use of the Heaviside step 
function 8: 



e{l+p~i) 



for I + p > i, 
for I + p < i. 



(C3) 



Similarly, the probability Dij of having i north coun- 
ters and j south counters with observed hits is given by 

l<i p^i 
m<j q<j 

Dij = X! X! ^IrnBpqflpifmqjQ {I + p - i) Q (m + q - j) , 

1=0 p=0 
m=0 q=0 

(C4) 

where Sim is the probability for having I north coun- 
ters and m south counters with signal hits and Bpq is 
the probability for having p north counters and q south 
counters with background hits. 



2. Derivation of Combinatoric Factor /; 



Ipi 



The combinatoric factor fipi gives the probability for 
observing i counters with hits given I counters with signal 
hits and p counters with background hits. The procedure 
described is based on two assumptions: 

1. Background hits are uncorrelated with signal hits. 
The background hits originate from interactions in 
earlier beam crossings, while the signal hits origi- 
nate from interactions in the current beam cross- 
ings, making this a reasonable assumption. 



2. The counters are either hit or not hit, ignoring the 
possibility that small amounts of charge from signal 
and background that are separately below thresh- 
old combine to exceed the discriminator threshold 
to form a hit. The same assumption is also made 
in the luminosity calculation and has been tested 
by comparing the north - south coincidence rate 
in a large sample of 10* simulated beam crossings 
with and without this assumption. No statistically 
significant difference was found. 

If u is the number of counters with both signal and 
background hits and v is the number of counters with 
only background hits, we obtain 



u = I + p — i, 
V = i — I. 



(C5) 



The number of ways to arrange u hits with both signal 
and background among I counters with signal hits is given 

by 



l\ 



u\{l-u)\ {l+p-i)\{i-p)V 



(C6) 



Similarly, the number of ways that v background-only 
hits can be distributed among the N — I counters without 
signal hits is 



N -i\ _ (N -ly. _ {N ~iy. 



(C7) 



Thus, the total number of background hit combinations 
is the product of Eqs. IC6l and lC7l 

l\fN-l\_ ^!(iV-0! 

uj\ V J ^ {i + p-iyii- py {i - ly {n ~ iy ■ 

(C8) 

The total number of ways to arrange p background hits 
among N counters is 



m 



N 



pj p\{N-py 



(C9) 



Each arrangement of background hits among the N coun- 
ters is assumed to be equally likely, therefore the proba- 
bility of having i counters with hits, given I counters with 
signal hits and p counters with background hits, is the 
number of arrangements meeting this condition divided 
by the total number of arrangements of the background 
hits: 



o 



fipi 



n{N^iy.p\ {N~py 



(CIO) 



{i+p- iy {i - py [i - ly. [n - iym 



While the derivation did not treat signal and background 
hits symmetrically, the final result includes the expected 
I O p symmetry. 
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3. Unfolding Procedure 

Once the data and background multiplicity distribu- 
tions are measured, an unfolding procedure is used to 
extract the signal multiplicity distribution. To simplify 
the notation, the 2D multiplicity distributions, with in- 
dexes ranging from to 24, are re-labeled to use a single 
index that ranges from to 624. For example, the 2D 
multiplicity distributions for counters with observed hits 
can be written as: 



Da = Dij where a = 25i + j 
Sp = Sim where /? = 251 + m 
B. 



(Cll) 



7 



Bpq where 7 = 25p - 



The tensor Tan-y is defined as follows: 



flpi fj 

.0 



mq 



for 



I <i,P<i,m<j,q< j, 



i<l+p,j<m + q 
otherwise. 



(C12) 

In this notation, Eq. IC4I for the convolution of signal 
and background multiplicity distributions can be written 
as 



Da — TafSjSjsB, 



7' 



(C13) 



where the convention of implied summation for repeated 
indices has been used. To solve for the signal multiplicity 
distribution in this notation, the folding and unfolding 
matrices F and U are defined: 



Fap 

U 



(C14) 
(C15) 



Substituting these matrices into Eq. IC13I the following 
matrix equations are extracted: 



D = FS 

S = P-^D = UD 



(C16) 
(C17) 



Thus, the signal multiplicity distribution is obtained 
by taking the product of the unfolding and data matrices. 
The inversion of the folding matrix F is only necessary 
for a simultaneous determination of the covariance ma- 
trix for the signal multiplicity distribution. Otherwise 
a linear equation solver can be used to find the signal 
multiplicity distribution S. 



Covariance Matrix for the Signal Multiplicity 
Distribution 



The unfolding procedure introduces correlations 
among the multiplicity bins, so that the uncertainty is 
represented by a 625 x 625 element covariance matrix. 
Standard error propagation techniques are used to deter- 
mine the covariance matrix for the unfolded signal mul- 
tiplicity distribution, 

i6D)^ = Tafiy (SS)^ B^ + Tap-^Sp (SB)^ . (C18) 
This result can be formulated as a matrix equation: 

Gay = TajSySp (^19) 

SD = F {6S) + G (SB) (C20) 
dS ^ F-^[{5D)~G{5B)] 

^U[iSD)~G{SB)]. (C21) 

The covariance matrix for the multiplicity distribution of 
counters with signal hits is given by: 

C^ = (^{SS) (SS) 

= (u [{SD) - G (SB)] [{SD) - G {SB)f 

(022) 

Since the observed signal -I- background multiplicity 
distribution is uncorrelated with the background multi- 
plicity distribution: 

= [/ ((SD) {5Df) + UG ({SB) (SBf) G^U'^ 



= UG^U^ + UGC''G' U' , 



TttT 



(C23) 



where and G^ are the covariance matrices for the 
data and background multiplicity distributions. The 
data and background covariance matrices simply contain 
diagonal binomial error terms: 

G^p - SapDa (1 - Da) /No (C24) 

Cf^ = Sat3Ba (1 - Ba) /Nb (C25) 

where 5ap is the Kronecker 5 function and Nd and Nb 
are the total number of entries in the data and back- 
ground multiplicity distributions, respectively. 



